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f~) I Abstract. The current status of big bang nucleosynthesis is summarized. Particular 

^ ■ attention is paid to recent observations of *He and ^Li and their systematic uncer- 

^^ ' tainties. Be and B are also discussed in connection to recent ^Li observations and the 

^— ,. , primordial ^Li abundance. 

,__! ! 1 Introduction 

>: 

lO i The simplicity of the standard model of Big Bang Nucleosynthesis (BBN) and 

r^ ■ its success when confronted with observations place the theory as one of the 

^T I cornerstones of modern cosmology. BBN is based on the inclusion of an extended 

^-^ ■ nuclear network into a homogeneous and isotropic cosmology. Apart from the 

^^ I input nuclear cross sections, the theory contains only a single parameter, namely 

^^ ■ the baryon-to-photon ratio, r). Other factors, such as the uncertainties in reaction 

fH I rates, and the neutron mean-hfe can be treated by standard statistical and Monte 

f) i' Carlo techniques |^,||J4||^. The theory then allows one to make predictions (with 

specified uncertainties) of the abundances of the light elements, D, ■^He, ^He, and 
''Li. As there exist several detailed reviews on BBN, I will briefly summarize the 
^ I key results and devote this contribution to the impact of recent observations of 

Cd ■ ^He and ^Li along with the related observations of Be and B. In referring to the 

^ I standard model, I will mean homogeneous nucleosynthesis, with three neutrino 

flavors {N^ = 3), and a neutron mean life of 886.7 ± 1.9 s |g]. 
, ^ The dominant product of big bang nucleosynthesis is *He, resulting in an 

jrt ■ abundance of close to 25% by mass. Lesser amounts of the other light elements 

are produced: D and ■^He at the level of about 10~^ by number, and ^Li at 
the level of 10"^" by number. The resulting abundances of the light elements 
are shown in Figure |l|, over the range in 7710 = IQ^'^j] between 1 and 10. The 
curves for the ^He mass fraction, Y, bracket the computed range based mainly 
on the uncertainty of the neutron mean-life. Uncertainties in the produced ^Li 
abundances have been adopted from the results in Hata et al.Q. Uncertainties 
in D and "^He production are small on the scale of this figure. The dark shaded 
boxes correspond to the observed abundances and will be discussed below. 

At present, there is a general concordance between the theoretical predictions 
and the observational data, particularly, for "^He and ^LiH. These two elements 
indicate that 77 lies in the range 1.55 < 77 < 4.45. There is limited agreement for 
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Fig. 1. The light element abundances from big bang nucleosynthesis as a function of 

mo- 

D/H as well, as will be discussed below. High D/H narrows the range to 1.5 < 
r] < 3.4 and low D/H is compatible at the 2cr level in the range 4.2 < rj < 5.3. 

2 Data 



2.1 



*He 



The primordial *He abundance is best determined from observations of Hell —> 
Hel recombination lines in extragalactic HII (ionized hydrogen) regions. There 
is a good collection of abundance information on the ^He mass fraction, Y, 0/H, 
and N/H in over 70 such regions [pp|,pi][. Since "^He is produced in stars along 
with heavier elements such as Oxygen, it is then expected that the primordial 



Big Bang Nucleosynthesis 



0.30 




200 



10^ O/H 

Fig. 2. The Hehum (Y) and Oxygen (O/H) abundances in extragalactic HII regions, 
from refs. |8| and from ref. [ho|. Lines connect the same regions observed by different 
groups. The regression shown leads to the primordial ''He abundance given in Eq. (hi). 



abundance of ''He can be determined from the intercept of the correlation be- 
tween Y and O/H, namely Yp = y(0/H ^ 0). A detailed analysis ||ll| of the 
data found 

Fp = 0.238 ± 0.002 ± 0.005 (1) 

The first uncertainty is purely statistical and the second uncertainty is an esti- 
mate of the systematic uncertainty in the primordial abundance determination. 
The solid box for ^He in Figure represents the range (at 2tTstat) from (nl). The 
dashed box extends this by including the systematic uncertainty. The He data 
is shown in Figure 0. 

The helium abundance used to derive dh) was determined using assumed 
electron densities n in the HII regions obtained from SII data. Izotov, Thuan, 
& Lipovetsky S proposed a method based on several He emission lines to "self- 
consistently" determine the electron density. Their data using this method yields 
a higher primordial value 



Yp = 0.244 ± 0.002 ± 0.005 



(2) 



As one can see, the resulting primordial ^He abundance shows significant sen- 
sitivity to the method of abundance determination, leading one to conclude that 
the systematic uncertainty (which is already dominant) may be underestimated. 
Indeed, the determination (|lj) of the primordial abundance above is based on a 
combination of the data in refs. g], which alone yield Yp — 0.228±0.005, and the 
data of ref. ||l^ (based on SII densities) which give 0.239 ±0.002. The abundance 
(||) is based solely on the self-consistent method yields and the data of |1^ . One 
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should also note that a recent determination |13] of the ''He abundance in a 
single object (the SMC) also using the self consistent method gives a primordial 
abundance of 0.234 ± 0.003 (actually, they observe Y = 0.240 ± 0.002 at [0/H] 
= -0.8, where [0/H] refers to the log of the Oxygen abundance relative to the 
solar value, in the units used in Figure 0, this corresponds to lO^O/H — 135). 
Therefore, it will useful to discuss some of the key sources of the uncertainties in 
the He abundance determinations and prospects for improvement. To this end, 
I will briefly discuss, the importance of reddening and underlying absorption in 
the H line line measurements, Monte Carlo methods for both H and He, and 
underlying absorption in He. 

The He abundance is typically quoted relative to H, e.g.. He line strengths 
are measured relative to H(3. The H data must first be corrected for underlying 
absorption and reddening. Beginning with an observed line flux F{X), and an 
equivalent width W{X), we can parameterize the correction for underlying stellar 
absorption as 

The parameter a is expected to be relatively insensitive to wavelength. A red- 
dening correction is applied to determine the intrinsic line intensity /(A) relative 
to HP 

where /(A) represents an assumed universal reddening law and C{Hj3) is the 
correction factor to be determined. By minimizing the differences between Xfi{X) 
to theoretical values, Xt(A), for A = Ha,H^ and H5, one can determine the 
parameters a and C{Hj3) self consistently [Q, and run a Monte Carlo over the 
input data to test the robustness of the solution and to determine the systematic 
uncertainty associated with these corrections. 

In Figure ^ |Q, the result of such a Monte-Carlo based on synthetic data 
with an assumed correction of 2 A for underlying absorption and a value for 
C{H(3) = 0.1 is shown. The synthetic data were assumed to have an intrinsic 
2% uncertainty. While the mean value of the Monte-Carlo results very accurately 
reproduces the input parameters, the spread in the values for a and C{Hf3) are 
considerably larger than one would have derived from the direct x^ minimization 
solution due to the covariance in a and C{H[3). 

The uncertainties found for H(3 must next be propagated into the analysis for 
^He. We can quantify the contribution to the overall He abundance uncertainty 
due to the reddening correction by propagating the error in eq. (0). Ignoring all 
other uncertainties in Xii{X) = I{X)/I{Hl3), we would write 

^^\nlOf{\)aciHm (5) 

In the example discussed above, (Jc(Hfi) ~ 0-04 (from the Monte Carlo), and 
values of / are 0.237, 0.208, 0.109, -0.225, -0.345, -0.396, for He lines at AA3889, 
4026, 4471, 5876, 6678, 7065, respectively. For the bluer lines, this correction 
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Fig. 3. A Monte Carlo determination of the underlying absorption a (in A ), and 
reddening parameter C{Hj3), based on synthetic data. 



alone is 1 - 2 % and must be added in quadrature to any other observational 
errors in Xji. For the redder lines, this uncertainty is 3 - 4 %. This represents the 
minimum uncertainty which must be included in the individual He I emission 
line strengths relative to H/3. 

Next one can perform an analogous procedure to that described above to 
determine the '^He abundance [|l3J . We again start with a set of observed quan- 
tities: line intensities /(A) which include the reddening correction previously 
determined along with its associated uncertainty which includes the uncertain- 
ties in C{H(3); the equivalent width W{X); and temperature t. The Helium line 
intensities are scaled to H(3 and the singly ionized helium abundance is given by 

+ /(A) EjHP) W{X)+aHei . 1 l_ f. 

y ^ ^ I{H[3) E{\) ^ W{\) \l+7)/(r) ^ ^ 

where E{\)/ E{H(3) is the theoretical emissivity scaled to H(3. The expression (||) 
also contains a correction factor for underlying stellar absorption, parameterized 
now by anei, a density dependent coUisional correction factor, (1 -1-7)"^, and a 
flourecence correction which depends on the optical depth r. Thus y^ implicitly 
depends on 3 unknowns, the electron density, n, anei, and r. 

One can use 3-6 lines to determine the weighted average helium abundance, 
y. From y, we can calculate the x^ deviation from the average, and minimize 
X^, to determine n, UHei, and r. Uncertainties in the output parameters are also 
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determined. In principle, under the assumption of small values for the optical 
depth t(3889), it is possible to use only the three bright lines A4471, A5876, and 
A6678 and still solve self-consistently for He/H, density, and anei- Of course, 
because these lines have relatively low sensitivities to collisional enhancement, 
the derived uncertainties in density will be large. 

The addition of A7065 was proposed Q as a density diagnostic and then, 
A3889 was later added to estimate the radiative transfer effects (since these 
are important for A7065). Thus the five line method has the potential of self- 
consistently determining the density and optical depth in the addition to the '*He 
abundance. The procedure described here differs somewhat from that proposed 
in ly , in that the x^ above is based on a straight weighted average, where as in [g[ 
the difference of a ratio of He abundances (to one wavelength, say A4471) to the 
theoretical ratio is minimized. When the reference line is particularly sensitive 
to a systematic effect such as underlying stellar absorption, this uncertainty 
propagates to all lines this way. 

Adding A4026 as a diagnostic line increases the leverage on detecting un- 
derlying stellar absorption. This is because the A4026 line is a relatively weak 
line. However, this also requires that the input spectrum is a very high quality 
one. A4026 is also provides exceptional leverage to underlying stellar absorption 
because it is a singlet line and therefore has very low sensitivity to collisional 
enhancement (i.e., n) and optical depth (i.e., t(3889)) effects. 

As in the case of the hydrogen lines, Monte-Carlo simulation of the He data 
can be used to test the robustness of the solution for n, anei, and r ||l3[. Figure 
Q presents the results of modeling of 6 synthetic He I line observations. The four 
panels show the results of a density — 100 cm~^, UHei = 0, and t(3889) = 
model. The solid lines show the input values (e.g., He/H = 0.080) for the original 
calculated spectrum. The solid circles (with error bars) show the results of the 
X^ minimization solution (with calculated errors) for the original synthetic input 
spectrum. The small points show the results of Monte Carlo realizations of the 
original input spectrum. The solid squares (with error bars) show the means and 
dispersions of the output values for the x^ minimization solutions of the A/Ionte 
Carlo realizations. 

Figure demonstrates several important points. First, the x^ minimization 
solution finds the correct input parameters with errors in He/H of about 1% 
(less than the 2% errors assumed on the input data, showing the power of using 
multiple lines). There is a systematic trend for the Monte Carlo realizations to 
tend toward higher values of Hc/H. This is because, the inclusion of errors has 
allowed minimizations which find lower values of the density and non-zero values 
of underlying absorption and optical depth. Note that the size of the error bars 
in He/H have expanded by roughly 50% as a result. We can conclude from this 
that simply adding additional lines or physical parameters in the minimization 
does not necessarily lead to the correct results. In order to use the minimization 
routines effectively, one must understand the role of the interdependencies of the 
individual lines on the different physical parameters. Here we have shown that 
trade-offs in underlying absorption and optical depth allow for good solutions at 
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Synthetic Data (T = 1.8, EW(HB) = 100) 
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Fig. 4. Results of modeling of 6 synthetic He I line observations. The four panels show 
the results of a density — 100 cni~^, uheI = 0, and r(3889) — model. 



densities which are too low and resulting in helium abundance determinations 
which are too high. Note that in the lower right panel of Figure Q that the values 
of the x^ do not correlate with the values of y+. The solutions at higher values 
of absorption and j/"*" are equally valid as those at lower absorption and y'^ . 

Figure ra shows the results of the Monte Carlo when both r and and 7^ 
0, and n — 100 cm""^. It is encouraging that in perhaps more realistic cases 
where the input parameters are non-zero, we are able to derive results very close 
to their correct values. The average of Monte Carlo realizations is remarkably 
close to the straight minimization for all of the derived parameters (n,aHei,T 
and y~^). However, there is an enormous dispersion in these results due to the 
degeneracy in the solutions with respect to the physical input parameters. This 
results in error estimates for parameters which are significantly larger than in the 
straight minimization. For example, the uncertainties in both the density and 
optical depth are almost a factor of 3 times larger in the Monte Carlo. When 
propagated into the uncertainty in the derived value for the He abundance, we 
find that the uncertainty in the Monte Carlo result (which we argue is a better, 
not merely more conservative, value) is a factor of 2.5 times the uncertainty 
obtained from a straight minimization using 6 line He lines. This amounts to 
an approximately 4% uncertainty in the He abundance, despite the fact that we 
assumed (in the synthetic data) 2% uncertainties in the input line strengths. This 
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Synthetic Data (T= 1.8, EW(HB) = 100) 
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Fig. 5. Similar plot to Figure W except that the underlying absorption is 0.1 A and 
t(3889) = 0.1. 



is an unavoidable consequence of the method - the Monte Carlo routine explores 
the degeneracies of the solutions and reveals the larger errors that should be 
associated with the solutions. 

In Figure |6[ I show the result of a single case based on the data of ref. jl^l for 
SBS1159+545. Here, the helium abundance and density solutions arc displayed. 
The vertical and horizontal lines show the position of the solution in ^^ . The 
circle shows the position of the our solution to the minimization, and the square 
shows the position of the mean of the Monte-Carlo distribution. The spread 



shown here is significantly greater than the uncertainty quoted in |10 



2.2 ^Li 

The abundance of ^Li has been determined by observations of over 100 hot, 
population-II halo stars, and is found to have a very nearly uniform abundance |14| . 
For stars with a surface temperature T > 5500 K and a metallicity less than 
about l/20th solar (so that effects such as stellar convection may not be impor- 
tant), the abundances show little or no dispersion beyond that which is consistent 
with the errors of individual measurements. The Li data from Ref. |15| indicate 
a mean ^Li abundance of 



Li/H= (1.6 ±0.1) X 10' 



10 



(7) 
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Fig. 6. A Monte Carlo determination of the helium abundance and electron density 
(in cm~^) for the region SBS11159+545. Solutions for a! and r are not shown here. 



The small error is statistical and is due to the large number of stars in which ^Li 
has been observed. The solid box for ^Li in Figure R^ represents the 2crstat range 
from (|). 

There is, however, an important source of systematic error due to the pos- 
sibility that Li has been depleted in these stars, though the lack of dispersion 
in the Li data limits the amount of depletion. In fact, a small observed slope 
in Li vs Fe and the tiny dispersion about that correlation indicates that deple- 
tion is negligible in these stars ||l^. Furthermore, the slope may indicate a lower 
abundance of Li than that in (6). The observation p7[ of the fragile isotope ^Li 
is another good indication that ^Li has not been destroyed in these stars |18[. 

The weighted mean of the ^Li abundance in the sample of ref. jl^] is [Li] 
= 2.12 ([Li] = log '^Li/H -\- 12) and is shghtly lower than that in eq. (0), the 
difference is a systematic effect due to analysis methods. It is common to test for 
the presence of a slope in the Li data by fitting a regression of the form [Li] = 
a + fi [Fe/H]. These data indicate a rather large slope, /3 = 0.07 — 0.16 and hence 

a downward shift in the "primordial" lithium abundance /\[Li] = —0.20 0.09. 

Models of galactic evolution which predict a small slope for [Li] vs. [Fe/H], can 
produce a value for /3 in the range 0.04 - 0.07 [|l9[. Of course, if we would like to 
extract the primordial ^Li abundance, we must examine the linear (rather than 
log) regressions. For Li/H = a' -I- ^'Fe/Fe©, we find a' = 1 — 1.2 x 10"^" and 
6' = 40 — 120 X 10^^". A similar result is found fitting Li vs O. Overall, when the 
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regression based on the data and other systematic effects are talcen into account 
a best value for Li/H was found to be |19[ 



Li/H = 1.23 X 10 



-10 



(8) 



with a plausible range between 0.9 - 1.9 xlO ^^. The dashed box in Figure Q 
corresponds to this range in Li/H. 
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Fig. 7. Contributions to the total predicted lithium abundance from the adopted GCE 
model of ref. pQ], compared with low metallicity stars (from pM) and a sample of high 
metallicity stars. The solid curve is the sum of all components. 

Figure shows the different Li components for a model with (^Li/H)p = 
1.23 X 10^^". The linear slope produced by the model is 6' = 65 x 10~^°, and is 
independent of the input primordial value (unlike the log slope given above). The 
model [ po| includes in addition to primordial ^Li, lithium produced in galactic 
cosmic ray nucleosynthesis (primarily a + a fusion), and ^Li produced by the 
i^-process during type II supernovae. As one can see, these processes are not suffi- 
cient to reproduce the population I abundance of ^Li, and additional production 
sources are needed. 



3 LiBeB 



The question that one should ask with regard to the above discussion of ^Li (and 
as we will see below when discussing concordance, ''Li will play an important 
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role in determining the baryon density 77), is whether or not there is additional 
evidence for the post big bang production of ^Li. There is in fact evidence in the 
related observation of the intermediate mass elements of ^Li, Be and B. While 
these elements are produced in the big bang [Q, their predicted primordial 
abundance is far below their observed abundance, which like ^Li is determined 
by observations of old metal poor halo stars. Whereas in the range 7710 = 1.5—4.5, 
standard BBN predicts abundances of 

^Li/H « (2 - 9) X 10"^'' 

9Be/H« (0.04-2) x 10"^^ 

ioB/H« (0.5-3) X 10-1^ 

"B/H « (0.02- 1) X 10"^^ (9) 

the observed abundances found in Pop II halo stars are: ^Li/H « few x 10~^^, 
^Be/H - 1 - 10 X 10^13^ and B/H ~ 1 - 10 x 10^^^ It is generally recognized 
that these isotopes are not of primordial origin, but rather have been produced 
in the Galaxy, through cosmic-ray nucleosynthesis. 

Be and B have been observed in the same pop II stars which show Li and 
in particular there are a dozen or so stars in which both Be and ^Li have been 
observed. Thus Be (and B though there is still a paucity of data) can be used as 
a consistency check on primordial Li. Based on the Be abundance found in these 
stars, one can conclude that no more than 10-20% of the ^Li is due to cosmic ray 
nucleosynthesis leaving the remainder (the abundance in Eq. (p|)) as primordial. 
This is consistent with the conclusion reached in Ref . |19|| . 

In principle, we can use the abundance information on the other LiBeB iso- 
topes to determine the abundance of the associated GCRN produced ^Li. As it 
turns out, the boron data is problematic for this purpose, as there is very likely 
an additional significant source for ^^B, namely i/-process nucleosynthesis in su- 
pernovae. Using the subset of the data for which Li and Be have been observed 
in the same stars, one can extract the primordial abundance of ^Li in the context 
of a given model of GCRN. For example, a specific GCRN model, predicts the 
ratio of Li/Be as a function of [Fe/H] . Under the (plausible) assumption that all 
of the observed Be is GCRN produced, the Li/Be ratio would yield the GCRN 
produced ^Li and could then be subtracted from each star to give a set of primor- 
dial ^Li abundances. This was done in [g2| where it was found that the plateau 
was indeed lowered by approximately 0.07 dex. However, it should be noted 
that this procedure is extremely model dependent. The predicted Li/Bc ratio 
in GCRN models was studied extensively in [g3|. It was found that Li/Be can 
vary between 10 and ^ 300 depending on the details of the cosmic-ray sources 
and propagation-e.g., source spectra shapes, escape pathlength magnitude and 
energy dependence, and kinematics. 

In contrast, the ^Li/^Li ratio is much better determined and far less model 
dependent since both are predominantly produced by a — a fusion rather than by 
spallation. The obvious problem however, is the paucity of ^Li data. As more ^Li 
data becomes available, it should be possible to obtain a better understanding 
of the relative contribution to ^Li from BBN and GCRN. 
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The associated BcB elements are clearly of importance in determining the 
primordial ^Li abundance, since Li is produced together with Be and B in accel- 
erated particle interactions such as cosmic ray spallation. However, these pro- 
duction processes are not yet fully understood. Standard cosmic-ray nucleosyn- 
thesis is dominated by interactions originating from accelerated protons and 
a's on CNO in the ISM, and predicts that BeB should be "secondary" versus 
the spallation targets, giving Be oc O^. However, this simple model was chal- 
lenged by the observations of BeB abundances in Pop II stars, and particularly 
the BeB trends versus metallicity. Measurements showed that both Be and B 
vary roughly linearly with Fe, a so-called "primary" scaling. If O and Fe are 
co-produced (i.e., if 0/Fe is constant at low metallicity) then the data clearly 
contradicts the canonical theory, i.e. BeB production via standard GCR's. 

There is growing evidence that the 0/Fe ratio is not constant at low metallicity | 
but rather increases towards low metallicity. This trend offers a solution to 
resolve discrepancy between the observed BeB abundances as a function of 
metallicity and the predicted secondary trend of GCR spallation [gO| . As noted 
above, standard GCR nucleosynthesis predicts Be ex O^, while observations show 
Be ~ Fe, roughly; these two trends can be consistent if 0/Fe is not constant in 
Pop II. A combination of standard GCR nucleosynthesis, and t^-process produc- 
tion of ^^B may be consistent with current data. 

Thus the nature of the production mechanism for BeB (primary vs. sec- 
ondary) rests with the determination of ratio of 0/Fe at low metallicity. In any 
case, it is clear that given a primary mechanism, it will be dominant in the early 
phases of the Galaxy, and secondary mechanisms will dominate in the latter 
stage of galactic evolution. The cross over or break point is uncertain. In Figure 
g, a plausible model for the evolution of BeB is shown and compared with the 
data Ig^ . Shown by the short dashed lines are standard galactic cosmic-ray nu- 
cleosynthesis, which is mostly secondary, but contains some primary production 
as well. The long dashed curves are purely primary, and in the case of boron, the 
u process has been included and this too is primary. The solid curves represent 
the total Be and B abundance as a function of [0/H]. As one can see such a 
model fits the data quite well. 



4 Concordance 

Let us now to turn to the question of concordance between the BBN predictions 
and the observations discussed above. This is best summarized in a comparison of 
likelihood functions as a function of the one free parameter of BBN, namely the 
baryon-to-photon ratio 77. By combining the theoretical predictions (and their 
uncertainties) with the observationally determined abundances discussed above, 
we can produce individual likelihood functions |^] which are shown in Figure ^ 
A range of primordial ^Li values are chosen based on the the abundances in Eqs. 
(0) and (g) as well as a higher and lower value. The double peaked nature of the 
Li likelihood functions is due to the presence of a minimum in the predicted 
lithium abundance in the expected range for rj. For a given observed value of 
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Fig. 8. Be vs O {top panel) and B vs O {bottom panel). Data shown are found to have 
a break point as indicated. Models are adjusted to achieve the break point and O/Fe 
slope of these data. 



^Li, there are two likely values of 77. As the lithium abundance is lowered, one 
tends toward the minimum of the BBN prediction, and the two peaks merge. 
Also shown are both values of the primordial ''He abundances discussed above. 
As one can see, at this level there is clearly concordance between ''He, ^Li and 
BBN. 
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Fig. 9. Likelihood distributions for four values of primordial ^Li/H (10^" x ^Li = 1.9 
(dashed), 1.6 (dotted), 1.23 (solid), and 0.9 (dash-dotted)), and for He (shaded) for 
which we adopt Yp = 0.238 ± 0.002 ± 0.005 (Eq. (1)). Also shown by the long dashed 
curve is the likelihood function based on the '*He abundance from Eq. (2). 



The combined likelihood, for fitting both elements simultaneously, is given 
by the product of two of the functions in Figure O. The combined likelihood is 
shown in Figure pG, for the two primordial values of ^Li in Eqs. m) and m). For 
^Lip = 1.6 X 10~ (shown as the dashed curve), the 95% CL region covers the 
range 1.55 < 7710 < 4.45, with the two peaks occurring at 7710 = 1-9 and 3.5. This 
range corresponds to values of Hb between 



0.006 < ilsh^ < .016 



(10) 



For ^Lip = 1.23 x 10" 



'^'^ (shown as the solid curve), the 95% CL region covers the 
range 1.75 < rjio < 3.90. In this case, the primordial value is low enough that 
the two lithium peaks are more or less merged as is the total likelihood function 
giving one broad peak centered at 7710 ~ 2.5. The corresponding values of fis in 
this case are between 



0.006 < ngh-^ < .014 



(11) 



When deuterium is folded into the mix, the situation becomes more compli- 
cated. Although there are several good measurements of deuterium in quasar 
absorption systems [|2^, and many of them giving a low value of D/H ~ (3.4 ± 
0.3) X 10^^ [^, there remains an observation with D/H nearly an order of 
magnitude higher D/H ~ (2.0 ± 0.5) x 10""* ^. 

Because there are no known astrophysical sites for the production of deu- 
terium, all observed D is assumed to be primordial. As a result, any firm deter- 
mination of a deuterium abundance establishes an upper bound on 77 which is 
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Fig. 10. Combined likeliiiood distributions for two values of primordial '^Li/H (10^° x 
■^Li = 1.6 (dashed), 1.23 (solid)), and ^He with Yp = 0.238 ± 0.002 ± 0.005 (Eq. (1)). 





1.6 




1.5 


^ 


1.4 


p- 


1.3 


3 


1 7 


°i 


1.1 


. 


1.0 


p 


0.9 


.-2 


0.8 


-"J 


0.7 




0.6 


. 


O.S 


p 


0.4 


„2 


0.3 


5 8 


0? 


"J 


0.1 



High D/H 




TllO 



Fig. 11. Likelihood distributions for two values of primordial ^Li/H (10^" x '^Li = 1.6 
(dashed) and 1.23 (solid)), and ''He with Yp = 0.238 ± 0.002 ± 0.005 from Eq. (Q) 
(shaded) and Yp = 0.244 ± 0.002 ± 0.005 from Eq. (||) (long dashed). Also shown are 
the two likelihood functions for high and low D/H as marked. 



robust. Thus the ISM measurements |2^ of D/H =1.6 xlO^^ imply an upper 
bound 7710 < 9. 

It is interesting to compare the results from the likelihood functions of ^He 
and "^Li with that of D/H. This comparison is shown in Figure |ll|. Using the 
higher value of D/H = (2.0 ± 0.5) x 10^^, we would find excellent agreement 
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between ''He, ^Li and D/H. The predicted range for 77 now becomes 



1.6 < 7710 < 3.2 



(12) 



with the peak likeHhood value at 7710 = 2.1, ''He and ^Li abundances from eqs. 



(m and (H) respectively. This corresponds to fish 



.008 



+ .004 
.002- 



The higher 



Li abundance of eq. (0) drops the peak value down slightly to ryio = 1.8 and 
broadens the range to 1.5 - 3.4. The higher ''He abundance shifts the peak and 
range (relative to eq. (p^)) up to 2.2 and 1.7 - 3.5. 
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Fig. 12. 50%, 68% & 95% C.L. contours of L47 and L247 where observed abundances 
are given by eqs. (nl and Vm, and high D/H. 



If instead, we assume that the low value of D/H = (3.4 ± 0.3) x 10 



IS 



the primordial abundance, there is hardly any overlap between the D and ^Li, 
particularly for the lower value of ^Li from eq. (^). There is also very limited 
overlap between D/H and "^He, though because of the flatness of the "^He abun- 
dance with respect to 77, as one can see, the likelihood function for the larger 
value of ''He from eq. (||) is very broad. In this case, D/H is just compatible 
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(at the 2 u level) with the other light elements, and the peak of the likelihood 
function occurs at roughly 7710 = 4.8 and with a range of 4.2 - 5.6. 
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Fig. 13. 50%, 68% & 95% C.L. contours of L47 and L247 where observed abundances 
are given by eqs. (nl and [7|), and low D/H. 



It is important to recall however, that the true uncertainty in the low D/H 
systems might be somewhat larger. Mesoturbulence efFects|Q allow D/H to be 
as large as 5 x 10~^. In this case, the peak of the D/H likelihood function shifts 
down to r/10 ~ 4, and there would be a near perfect overlap with the high r\ 
''Li peak and since the ^He distribution function is very broad, this would be a 
highly compatible solution. 

We can obtain still more information regarding the compatibility of the 
observed abundance and BBN by considering generalized likelihood functions 
where we allow N^, to vary as well |3,^,^,Q. The likelihood functions now 
become functions of two parameters £(77, Nv)- 

The peaks of the distribution as well as the allowed ranges of r\ and 7V^ are 
easily discerned in the contour plots of Figures |l2| and |3| which show the 50%, 
68% and 95% confidence level contours in L47 and L247 projected onto the 77- 



18 Keith A. Olive 

N^ plane, for high and low D/H as indicated. L47 corresponds to the likelihood 
function based on ''He and ^Li only, whereas L247 includes D/H as well. The 
crosses show the location of the peaks of the likelihood functions. L47 peaks at 
N^ = 3.2, 7710 = 1.85 and at N^ = 2.6, 7710 = 3.6. The 95% confidence level 
allows the following ranges in rj and N^, 

1.7 <N^< 4.3 1.4 < ?7io < 4.9 (13) 

Note however that the ranges in ry and N,j are strongly correlated as is evident 
in Figure 02. 

With high D/H, ^247 peaks at Ny — 3.3, and also at r/10 = 1.85. In this case 
the 95% contour gives the ranges 

2.2 <Ny< 4.4 1.4 < r/10 < 2.4 (14) 

Note that within the 95% CL range, there is also a small area with 7710 = 3.2 — 3.5 
and Ny = 2.5 - 2.9. 

Similarly, for low D/H, L247 peaks at N^ = 2.4, and 7710 = 4.55. The 95% 
CL upper limit is now iVj^ < 3.2, and the range for 77 is 3.9 < 7710 < 5.4. It is 
important to stress that these abundances are now consistent with the standard 
model value of N^ = 3 at the 2 a level. 
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